LIEB-ROBINSON BOUNDS IN QUANTUM MANY-BODY PHYSICS 



BRUNO NACHTERGAELE AND ROBERT SIMS 



Abstract. We give an overview of recent results on Lieb-Robinson bounds and some of 
their applications in the study of quantum many-body models in condensed matter physics. 



1. Introduction 

Condensed matter physics and quantum computation have a common interest in 'comph- 
cated' states of quantum many-body systems. Low-temperature physics is often described 
by states with strong correlations and many interesting physical properties are due to those 
correlations. Quite often, these correlations are intimately related to fundamental features 
of quantum mechanics in a phenomenon that is usually referred to as entanglement. In 
the case of quantum computation entanglement is the essential feature without which the 
subject would not exist. 

The aim of these lecture notes is to present recent results that emerged from this shared 
interest in entangled quantum states. Entanglement is often described as 'non-local corre- 
lations' that are of a fundamentally different nature than correlations in the classical sense. 
The latter can be described as a property of a multivariate probability distribution. While 
it is absolutely correct that quantum entanglement is something new that is not found in 
any classical description of physics, it would not be correct to say that the complexity im- 
plied by entanglement precludes a useful discussion of locality in quantum physics. One 
of the key tools to unravel the complexity of typical ground states of a condensed matter 
physics model, states that could potentially be of interest for the implementation of a quan- 
tum computer, are the so-called Lieb-Robinson bounds. Lieb-Robinson bounds imply that 
non-relativistic quantum dynamics has, at least approximately, the same kind of locality 
structure provided in a field theory by the finiteness of the speed of light. 

The original work by Lieb and Robinson dates back to 1972 [29]. Since the work of 
Hastings [21j there have been a series of extensions and improvements [38 1 [23 | [35l [36l HO l H2] . 
These extension were largely motivated by the possibility of applying Lieb-Robinson bounds 
to a variety of problems concerning correlated and entangled states of quantum many-body 
systems. Other applications are focused on the dynamics of such systems. For example, a 
Lieb-Robinson type estimate was used in [12] to derive the nonlinear Hartree equation from 
mean-field many body dynamics. In [37], the application is to prove the existence of the 
dynamics of a class of anharmonic oscillator lattices in the thermodynamic limit. 

In these lecture notes we will discuss Lieb-Robinson bounds, also called locality or propa- 
gation estimates, depending on one's point of view, for quantum lattice systems and some of 
the applications in which they play a central role. Section [2] is devoted to the Lieb-Robinson 
bounds themselves. We start by considering systems with bounded interactions, of which 
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quantum spin systems are the typical examples, but which also includes models with an 
infinite-dimensional Hilbert space on each site, as long as the only unbounded terms in the 
Hamiltonian act on single sites. An example of this situation is the quantum rotor model. 
Next, in Section [2. 2 1 we study a class of systems with unbounded interactions: the harmonic 
oscillator lattice models and suitable anharmonic perturbations of them. 

In Section [3] we show how Lieb- Robinson bounds can be used to prove the existence of the 
dynamics in the thermodynamic limit. Intuitively, it is easy to understand that a locality 
property of some kind is essential for the dynamics of an infinite system of interacting 
particles to be well-defined. Again, we first treat systems with bounded interactions only, 
and then turn to the anharmonic oscillator lattices systems for which new subtleties arise. 

The first direct application to ground states of quantum many-body systems is the expo- 
nential clustering theorem which, succinctly stated, says that a non-vanishing spectral gap 
above the ground state implies the existence of a finite correlation length in that ground 
state. We sketch a proof of this results in Section HI One can obtain much more detailed 
information about the structure of the ground state. A good example of this is the Area 
Law for the entanglement entropy, which Hastings proved for one-dimensional systems [22] . 
We have made some progress to understand the structure of gapped ground states in higher 
dimensions, but a general Area Law is still lacking. The structure of gapped ground states 
that is now revealing itself can be regarded as a generalization of a special type of states that 
were shown to be the exact ground states of particular Hamiltonians, called valence bond 
solid models, of which the first non-trivial examples were introduced by Affleck, Kennedy, 
Lieb, and Tasaki in 1987 [2], most notable the spin-1 chain now called the AKLT model. 
These examples were generalized considerably in [14], where the class of finitely correlated 
states was introduced, a special subclass of which was later dubbed matrix product states. 
Since it is now becoming clear that the special structure of these states is approximately 
present in all gapped ground states, we provide in this section also a brief overview of the 
AKLT ground state. 

New applications and extensions of Lieb-Robinson bounds and techniques that employ 
them continue to be found. In Section [5l we briefly mention a few examples of recent results 
that are not discussed in detail in these notes. 

2. LiEB-RoBiNSON Bounds 

2.1. Bounded Interactions. In this section, we will consider general quantum systems 
defined on a finite sets T equipped with a metric d. We introduce these systems as follows. 

To each site x G T, we will associate a Hilbert space Tlx- In the context of quantum spins 
systems, the Hilbert space Ux is finite dimensional, whereas for oscillator systems, typically 
T-Lx = L'^i^^dqx)- For the locality results we present in this section, both systems can be 
treated within the same framework. 

With any subset X CT, the Hilbert space of states over X is given by 

(2.1) -Hx = (g)^x-, 

and the algebra of local observables over X is then defined to be 

(2.2) = (g) BiTix), 

where B{Hx) denotes the algebra of bounded linear operators on Hx- 
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The locality results we prove are expressed in terms of the support of local observables. 
Here, the support of an observable is understood as follows, li X CT, wc identify A G Ax 
with A (gi 1y\x £ -4r- In a similar manner, we have that for each X C Y G T, Ax C Ay- 
The support of an observable A G .4r is then the minimal set X C F for which A = A'(g)]lp\x 
with A' e Ax- 

The models we consider here correspond to bounded perturbations of local self-adjoint 
Hamiltonians. Specifically, we will fix a collection of on-site local operators {Hx}xer where 
each Hx is a sclf-adjoint operator with domain D{Hx) C T-Lx- In addition, wc will consider 
a general class of bounded perturbations. These are defined in terms of an interaction $, 
which is a map from the set of subsets of T to with the property that for each set X (ZV, 
$(X) G Ax and = 

Lieb-Robinson bounds are essentially an upper bound for the velocity with which per- 
turbations can propagate through the system. The metric d measure distances in the 
underlying space and it turns out that some regularity of V has to be assumed, which can 
be interpreted as a condition which guarantees that F can be 'nicely' embedded in a finite- 
dimensional space. This property is expressed in terms of a non-increasing, real-valued 
function F : [0, oo) — > (0, oo), which will enter in our estimate of the Lieb-Robinson ve- 
locity. F will also be used to impose a decay condition on the interactions in the system. 
The existence of a function F with the required properties is non-trivial only when the 
cardinality of F is infinite, but all the relevant quantities can be defined for a finite system. 
To each pair x, y G F, there is a number Cx,y such that 

(2.3) ^ (^(^' ^ ('^(^' y)) ^ ^-'V^ ('^(^' y)) ■ 

Take Cx,y to be the infimum over all such Cx,y, and denote by C = mB;Kx,y&T Cx,y Explicitly, 
for models with F C Z*^, one choice of F is given by F{r) = (1 -|- rY+'^. In this case, the 
convolution constant may be taken as C = 2*^+^ Z^^er -^(I^D- -'-^ general, the quantity 

^ ||<I>(A-)[| 
;yev ^ F{d{x,y))' 
x,yeX 

which is finite for any interaction ^> over F, will also play a role in our analysis. 

Now, for a fixed sequence of local Hamiltonians {Hx}xeT^ as described above, an inter- 
action and any subset A C F, we will consider sclf-adjoint Hamiltonians of the form 

(2.5) i^A = /fi- + Ht = Y.H.+ Y. ^ 

with domain ^xeA ^(^^) ^ Since each operator i^A is self-adjoint, it generates 

a Heisenberg dynamics, or time evolution, {r^^}, which is the one parameter group of 
automorphisms defined by 

(2.6) Tt^{A) = e^tHA j^^-itHj, fQj, g^^y ^ g 

Let us consider two concrete models of the type described above. 



(2.4) 
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Example 2.1. Take T-Lx = for each x G F and consider a Heisenberg Hamiltonian of 
the form 

(2.7) Hr = ^ BSx + ^ JxySx ■ Sy . 

d{x,y)<l 

As we will see below, the velocity corresponding to such a model depends on the interaction 
strength Jxy, and it is independent of the external magnetic field B. 

Example 2.2. Take T-Lx = L^(M, dqx) for each x E F and consider an oscillator Hamiltonian 
of the form 

(2.8) Hr = Y,Pl + y E ^^^--%)- 

d{x,y)<l 

If Hx = Px + y{Qx) is self-adjoint, then Theorem \2.3\ below estimates the velocity of such a 
model for any real-valued <j) £ L°°(IR). 

For a Hamiltonian that is a sum of local interaction terms, nearest neighbor for example, 
it is reasonable to expect that the spread of the support of a time evolved observable depends 
only on the surface area of the support of the observable being evolved; not the full volume 
of the support. This is important in some applications, e.g., to prove the split property 
of gapped spin chains [33J. To express the dependence on the surface area we will use the 
following notation. Let X C A C F. Denote the surface of X, regarded as a subset of A, by 

(2.9) Sa{X) = {Z cA : Z n X / and Z n 0} , 

and we will write S{X) = Sr{X). The <I>-boundary of a set X, written d^X, is given by 

(2.10) d^X = {x € X : 3Z € S{X) with x G Z and ^{Z) / } . 

Denote the distance between two sets AT, y C F by d{X,Y) = inmx£x,yeY d{x,y). 

The main result of this section is the following theorem, which was first proved in [36] . 

Theorem 2.3. Let T be a finite set and fix a collection of local Hamiltonians {Hx}x&r o,nd 
an interaction $ over F. Let X and Y be subsets ofV with d{X,Y) > and take any set 
A D X UY . For any pair of local observables A S Ax and B G Ay, the estimate 

(2.11) \\[rHA\B]\\ < ?MM (e^cil^lll^l _ j^^j^^y), 
holds for all t G M. Here 

(2.12) D{X,Y) = mm 



^ 5]F(d(x,y)),^ F{d{x,y)) 



Before we prove Theorem 12.31 we make a few comments which may be useful in inter- 
preting this result. First, we note that if X and Y have a non-empty intersection, then the 
argument provided below produces an analogous bound with the factor 11*11 *^ 1*1 — 1 replaced 
by 11*11 l*L In the case of empty intersection and for small values of \t\, (|2.1ip is a better 
and sometimes more useful estimates than the obvious bound || [rt(yl), 2?] || < 2||j4||||i?||, 
valid for all t G M. 

Next, (12. lip provides a fairly explicit locality estimate for the corresponding dynamics. 
We have thus far expressed this in terms of the function F defined above. Note that if F 
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is equipped with such a function F, then for every ^ > 0, F^{r) = e ^'^F{r) also satisfies 
(|2.3p above. Setting then 

(2.3, ll*t^,;-E 



^■'S/er Ff,{d{x,y)) 

x,yeX 



for all ;U > 0, it is easy to see that 



(2.14) Df,{X,Y) < ||$||^ min(|a$X| ,|a$y|) e-^'^^^'^) maxV F((i(x,?/)). 

In this case, ()2.1ip implies 
(2.15) 

||[r/^(^),5]|| < ^**^Jj**^** maxVF(d(x,y)) mm{\d^X\ ,\d^Y\) ^-^'['i(^^Y)-'-M^\t\] ^ 

i.e. the locality bounds decay exponentially in space with arbitrary rate > 0. For every 
/i, the system's velocity of propagation, satisfies the bound 



(2.16) u$ < 

A* 

As a final comment, we observe that for fixed local observables A and B, the bounds 
above, ()2.1ip and ()2.15p . are independent of the volume A C T; given that A contains the 
supports of both A and B. Furthermore, we note that these bounds place a constraint 
on only the minimum of the support of the two observables. Thus the estimate is still 
independent of A even if the support of one of the observables depends on the volume A. 

The proof of Theorem 12.31 uses a basic lemma about the growth of the solutions of first 
order, inhomogeneous differential equations. We state and prove it before the proof of the 
theorem. 

Let B he a, Banach space. For each f G M, let A{t) : B ^ B he a linear operator, and 
denote by X{t) the solution of the differential equation 

(2.17) dtX{t) = A{t)X{t) 

with boundary condition ^(0) = xq £ B. We say that the family of operators A(t) is 
norm-preserving if for every xq £ B, the mapping ■ B ^ B which associates xq i— t- X(t), 
i.e., 7t(a;o) = X{t), satisfies 

(2.18) II 7t(xo) II = II Xq II for all t e M. 

Some obvious examples are the case where ;B is a Hilbert space and A(t) is anti-hermitian 
for each t, or when ;B is a *-algebra of operators on a Hilbert space with a spectral norm 
and, for each t, A{t) is a derivation commuting with the *-operation. 

Lemma 2.4. Let A{t), for t € M, be a family of norm preserving operators in some Banach 
space B. For any function i? : M — )• 5, the solution of 

(2.19) dtY{t) = Ait)Y{t) + Bit), 
with boundary condition Y{0) = yQ, satisfies the bound 

(2.20) \\Y{t) - 7t(yo)|| < f\\B{t')\\dt'. 

Jo 
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Proof. For any t £ M, let X{t) be the solution of 

(2.21) dtX{t) = A{t)X{t) 

with boundary condition -'^(O) = xq, and let jt be the linear mapping which takes xq to 
X{t). By variation of constants, the solution of the inhomogeneous equation (|2.19p may be 
expressed as 

(2.22) Y{t) = 7t (yo + l^^s)-' {B{s))ds^ . 

The estimate (12.201) follows from (I2.22P as A{t) is norm preserving. □ 

Proof of Theorem \2.3[ Fix A C F as in the statement of the theorem. As this set will remain 
fixed throughout the argument, we will suppress it in our notation. In particular, we will 
denote r^^ merely by rj. 

Without loss of generality, we will assume that 

(2.23) D{X,Y)= Y^(d{x,y)). 

x£d^X yeY 

Otherwise, we apply the argument below to ||[T_f(i?),^]||. 
For each Z C F, we introduce the quantity 

(2^24) C,(Z:t):= sup i^M. 

and note that Cij(Z;0) < 2\\B\\6y{Z), where we defined 6y{Z) = 1 if y n Z 7^ and 
Sy{Z) = if y n Z = 0. Note that the dynamics generated by 

zeA zcx 
leaves Ax invariant. More precisely, if we define tI°^ by 

(2.25) ri-(^) = e^«^+^*)^e-<^A^+^*) for ah ^ G ^a, 

then we have that for every A G Ax, Tt°'^{A) G Ax for all t G M. This implies, recalling the 
definition ^JMi that 

(2.26) CB{X;t) = sup r--- . 

AeAxA^O \\A\\ 

Consider the function 

(2.27) fit) := [T,(ri-(.4)),i?], 

for A G Ax, B G Ay, and t G M. It is straightforward to verify that 

(2.28) f'{t) = ^ [rtmZ))J{t)]-i ^ [nir'^f (A)) , [n (HZ)) , B] . 

Zg5a(X) ZGSAiX) 

The first term in the above differential equation is norm preserving, see Lemma 12. 4| and 
therefore we have the bound 

(2.29) ||/(t)|| < 11/(0)11 + 2PII Y r\\[rs{HZ)),B]\\ds. 

zes(x) •'^ 
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Recalling definition ()2.24p . the above inequality readily implies that 

r\t\ 

(2.30) CB{X,t)<CB{X,0)+2 mZ)\\ CB{Z,s)ds, 

ZeS{X) •'^ 

where we have used (|2.26p . Iteration of (|2.30p yields that 

(2.31) CB{X,t) < 2\\B\\ ^^i-l^an, 



n! 

71=1 



where for n > 1, 



(2.32) an= E E E 5y{Zn)\{mZ,)\\ 

Zi€SiX) Z2eS(Zi) z„es(z„-i) i=i 

For any interaction <I>, one may estimate that 

(2.33) «i < E E 11^(^)11 ^ ii'^iiE E pidi^^y))- 

In addition, 



y&Y Z!=S{X): yeY xGdiSfX 

yez 



a2 < H E E 11^(^2)11 

y£Y Zi&S{X) zi£di,Zi z^cr: 

zi,y£Z2 

ziGZi 

< ll'^'ll' E YJ2P^d{x,z,))F{d{zuy)) 

x£d^X yeY 2iGr 

(2.34) < m^C Y Y^(d{x,y)), 

xed^X y&Y 

where we have used C from (|2.3p for the final inequality. With analogous arguments, one 
finds that for all n > 1, 

(2.35) an<mrC^'-' Y Y^^dix,y)). 

x&d<s>X y&Y 

Inserting (|2.35p into ()2.3ip we see that 

(2.36) CBiX,t) < ^ (e^^H*lll*l - l) ^ E P(dix,y)) , 

xed^X yeY 

from which ()2.1ip immediately follows. □ 

2.2. Unbounded Interactions. In this section we will consider locality estimates for sys- 
tems with unbounded interaction terms. There are very few results in this context. Results 
bounding the speed of propagation of perturbations in classical anharmonic lattice systems 
have been obtained [32l EH [H] , but these works do not provide explicit estimates for the 
Lieb-Robinson velocity. For a class of classical models similar to the quantum models we 
will discuss here, bounds for the Lieb-Robinson velocity have been proved recently in [42] . 
The only known results for quantum systems in this context apply to harmonic systems and 
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a class of bounded perturbations and were obtained in [36] and |37j . These are the results 
we will review here. 

We begin in Section 12.2.11 by introducing a well-known family of harmonic oscillator 
models defined on finite subsets of Z"'; compare with Example 12.21 Then, in Section [2.2.2^ 
we introduce a convenient class of observables, the Weyl operators, which the harmonic 
dynamics leaves invariant. In Section 12.2.31 we demonstrate a Lieb- Robinson bound for 
these harmonic models. Finally, in Section 12.2.41 we show that a similar Lieb-Robinson 
bound holds for a large class of bounded perturbations. 

2.2.1. Harmonic Oscillators. We first consider a system of coupled harmonic oscillators 
restricted to a finite volume. Specifically on cubic subsets A^, = {—L, L]^ C Z'^, we analyze 
Hamiltonians of the form 

d 

(2.37) Hi= Y,pI + ^\1 + Y. - 9-+-.)' 
acting in the Hilbert space 

(2.38) -Ha^ = (g) L\R,dq,). 

xeAL 

Here the quantities px and Qx, which appear in (j2.37|) above, are the single site momentum 
and position operators regarded as operators on the full Hilbert space TiA^ by setting 

(2.39) = 11(g) ••• (g) If^ -i-^ (g) 1- •• (g) 1 and = 1 (g • • • (g 1 (g g (g 1 • • • (g 1, 

dq 

i.e., these operators act non-trivially only in the x-th factor oH-Lal- These operators satisfy 
the canonical commutation relations 

(2.40) \Px,Py\ = [qx,qy\ = and [qx.Py] = iSx,y, 

valid for all x,y G A^. In addition, {ej}j^i are the canonical basis vectors in Z"^, the 
numbers Xj > and uj > are the parameters of the system, and the Hamiltonian is 
assumed to have periodic boundary conditions, in the sense that qx+e^ = Qx-{2L-i)ej if 
X £ Al but X + Cj ^ Al- It is well-known that these Hamiltonians are essentially self- 
adjoint on Cq°, see e.g |43) . Moreover, these operators have a diagonal representation in 
Fourier space. We review this quickly to establish some notation and refer the interested 
reader to [36j for more details. 
Introduce the operators 

(2.41) Qk = ^= ^ e-^'-^qx and = — ^ J] e'^^-p, , 
defined for each k £ A*^ = : xgAl}, and set 

(2.42) 7(A;) 



d 

a;2 + 4^Aj sin^{kj/2). 



A calculation shows that 

(2.43) = E ^(^) (2^fe&fe + 1; 
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where the operators bk and 6^ satisfy 



(2.44) bk = / --Pk-i\I^Qk and bl = -. P-k + Q-k ■ 

In this sense, we regard the Hamiltonian as diagonahzable. The special case of w = 
is discussed in [361. 



2.2.2. Weyl Operators. For these harmonic Hamiltonians, a specific class of observables, 
namely the Weyl operators, is particularly convenient. Given any function / : A^^ — t- C, the 
corresponding Weyl operator W{f) is defined by setting 



(2.45) W{f) = exp 



i ^ Ke[f{x)]qx + lva[f{x)]px 



It is easy to see that each W{f) is a unitary operator with 

(2.46) Wifr' = W{-f) = W{fr. 
Moreover, using the well-known Baker-Campbell-Hausdorff formula 

(2.47) e^+^ = e^e^e-[^-^]/2 if [.4, [A, B]] = [B, [A, B]] = , 

and the commutation relations ()2.40l) . it follows that Weyl operators satisfy the Weyl rela- 
tions 

(2.48) W{f)W{g) = W{g)W{f)e~'^'^^^f'9)] = w{f + g)e-^^mf,9)]/2 ^ 

for any /, (7 : A^, — )■ C. These operators also generate shifts of the position and momentum 
operators in the sense that 

(2.49) VF(/) = - Im[/(x)] and W{frp^W{f) = + Re[f{x)] . 

The algebra generated by all such Weyl operators is called the Weyl algebra. 

A key observation which we will exploit in our locality estimates is the fact that the 
harmonic dynamics leaves the Weyl algebra invariant. We state this as a lemma. Fix 
L > 1 and t S M. Denote by t^'^ the harmonic dynamics generated by H^, i.e. for any 
A e BiTiA^) set 

(2.50) rf'^(A) = e^*^Ue-^*^£ . 

Lemma 2.5. Fix L > 1 and t £ M.. There exists a mapping T^'^ : ^^(A^) — )• £^(A/^) for 
which 

(2.51) n'\wif)) = w^(r^^/) 

for any f £ i'^iAi). 

Proof. Since we will fix L > 1 throughout this argument and only consider harmonic Hamil- 
tonians, we will denote by = t^'^ and Tt = T^'^ to ease notation. To prove this lemma, 
it is also convenient to express a given Weyl operator in terms of annihilation and creation 
operators, i.e., 

(2.52) = {qx + iPx) and a* = {q^ - ip^) , 
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which satisfy 

(2.53) [a^, Uy] = [a*, a*] = and K., a*] = 5x,y for ah x,y £ Al- 

One finds that 



(2.54) W{f) = exp 

where, for each / G £'^{Al), we have set 
(2.55) 



— iaif) + a*{f)) 



«(/) = E ■/■(^)"- = E /(^x- 

xgAi i^eAj, 

It is easy to see that the harmonic dynamics acts trivially on the diagonalizing operators 
b, i.e., 

(2.56) rtibk) = e-^'-'^'^Hk and tM) = e^'^'^'Hl , 

where bk and 6^ are as defined in (j2.44p . Hence, if we further introduce 



(2.57) 



for each x € Al and, analogously to (I2.55p . define 

(2.58) b{f) = Y.W)b,, b*{f) = f{x)b:, 

for each / G £^(Al), then one has that 

(2.59) Tt{b{f)) = b{[T-HltT]f) , 

where J- is the unitary Fourier transform from £^(Ai) to ^^(A2)and Mt is the operator of 
multiplication by e^*'''^^^* in Fourier space with ^{k) as in (j2.42p . The proof is completed by 
demonstrating a change of variables relation between the a's and the Vs. 

A short calculation shows that there exists a linear mapping U : ^^(A^,) — )■ ^^(A^,) and 
an anti- linear mapping V : £'^{Al) — ). £'^{Al) for which 

(2.60) b{f)=aiUf) + a*iVf), 

sl relation know in the literature as a Bogoliubov transformation [30]. In fact, one has that 



(2.61) 



U = -F-^My,F and V = -F-'Mt_TJ 



-1 



2 + 2 

where J is complex conjugation and Mr^ is the operator of multiplication by 

1 



(2.62) 



r±(fc) 



again, with 7(/c) as in ()2.42p . Using the fact that F-i- is real valued and even, it is easy to 
check that 



(2.63) 

and 

(2.64) 



U*U - V*V = 1 = uu* - vv* 



V*U - U*V = = vu* - uv* 
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where we stress that V* is the adjoint of the anti-linear mapping V. The relation (j2.60p is 
invertible, in fact, 



(2.65) 

and therefore 
(2.66) 

Clearly then, 
(2.67) 



a{f) = b{U*f)-b*{V*f), 



Wif) = exp 



V2 



{b{iU* - V*)f) + b*{{U* - V*)f)) 



TtiWif)) = W{Ttf) , 
where the mapping Tt is given by 

(2.68) Tt = {U + V)F-^MtF{U* - V*) , 

and we have used (j2.59p . 



□ 



2.2.3. Lieb-Robinson bounds for harmonic Hamiltonians. In this section, we demonstrate a 
Lieb-Robinson type bound for these harmonic lattice systems. We state our estimate, first 
proved in [36|, as follows. 



Theorem 2.6. Fix L > 1. For any n > 0, the estimate 

(2.69) \\[r^^'^{W{f)),W{g)]\\<C J] \f{x)\\g{y)\e 

holds for all functions f,g G £^{Al) and any t G M . Here 

d 

(2.70) d{x, y) = / min \xj — yj + 2Lr]j \ . 

is the distance on the torus. Moreover 

(2.71) C = (l + c^,,e^/' + cj_\) 
mt/iCa,,A = (u;2 + 4EJ=iA,)V2. 

Before we prove this result, we make a few comments. First, we denote by 



(2.72) 



^;/,(^) = c^,Amax(^^,e(^/2)+i^ 



our estimate on the harmonic Lieb-Robinson velocity corresponding to decay rate ^ > 0. 
Optimizing over /u > produces a rate 1/2 < /io < 1 for which Vhino) < 4ctj^A- Next, in [36], 
see also [37], it is shown that the mapping t/*'^ appearing in Lemma 12.51 can be expressed 
as a convolution. In fact, dropping the dependency on h and L, one has that 



(2.73) 



Ttf = f*{h\ 



,(0) 



-Ahr>+h^^)]+jj'-ih^^-ht'^: 
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where 



(2.74) 



h[ ^\x) = Im 



h["\x) = Re 



h[^\x) = Im 



J_ 



J,k-x—2i'y(k)t 



1 ik-x-2ij{k)t 



feeAt 



jife-x— 2i7(fc)t 



|Al| 



fceA* 



Next, by direct calculation, the following is proven in [36]. 

Lemma 2.7. Consider the functions defined in \2. 74 )■ For w > 0, Ai, . . . , A,^ > 0, but such 
that c^^x = (cj^ + ^jY^"^ > 0, and any fJ- > 0, the bounds 



(2.75) 



hfix) 
h?{x) 



1 -M(|x|-c„,Amax(|,eW2)+i)|i|) 



< c^,xe>^'^e 



/./2 -m(|^|-c.,a max(| , e(''/2)+l) \t\) 



hold for all t and x G A^,. Here \x\ = Yl'j=i 

Proof of Theorem \2.(k With L > 1 fixed, we again drop the dependence of the dynamics 
on h and L. Observe that 

[Tt{W{f)),W{g)] = {W{Ttf)-W{g)W{Ttf)W{-g)}W{g) 

(2.76) = {l-e*^™[<^*-^'^'>l}H^(rt/)Ty(5), 

where we have used Lemma 12.51 and the Weyl relations (I2.48P . Since Weyl operators are 
unitary, the estimate 

(2.77) \\[n{W{f)),W{g)]\\ < \lm[{Ttf,g)]\ 
readily follows. The bound 

(2.78) \lm[{Ttf,g)]\ < \Ttf{y)\\g{y)\ 



S/GAi 



□ 



is obvious and (|2.69|) is now a consequence of (|2.73|) and Lemma [221 

In analogy with Section [2. !( these Lieb- Robinson bounds can be expressed in terms of a 
family of non-increasing, real- valued functions : [0, oo) — )• (0, oo), parametrized by // > 0, 
given by 

(2-79) W=(YT7)^- 

With d{-,-) the metric on the torus, see (I2.83p . it is easy to see that 

(2.80) Yl F^{d{x,z))F^{d{z,y)) < CdF^{d{x,y)) 

zeAL 
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with 

(2-81) Cd = 2^^+^ V 7 ^--r-^ . 



The fohowing corollary of Theorem 12.61 is immediate. 

Corollary 2.8. Fix L > 1. For any /i > and e > 0, the estimate 

(2.82) ||[r,'''^t^(/)),W(5)]|| < C(e,^)e('^+^>'>('^+^)l*l ^ \f{x)\\g{y)\ F,{d{x,y)) 

holds for all functions f,g£ £'^{Al) and any t gW. Here 

d 

(2.83) d{x, y) = y min \xj — yj + 2Lr]j \ . 
is the distance on the torus, 

(2.84) C(e,M) = (l + c^,Ae(^+^)/' + cj\) supe-"(l + s)^+' , 

Cu},x = (f^^ + 4 l^}'=i -^j)^^^; o.nd VhifJ-) = Ci_j^x max(^, e^/^"'"^) is the harmonic velocity corre- 
sponding to the decay rate /j,. 

2.2.4. Lieb-Robinson bounds for anharmonic systems. In this section we will consider per- 
turbations of the harmonic Hamiltonians defined in Section 12.2.11 The results we prove here 
originally appeared in [36j and [37]. The perturbations are defined as follows. Fix L > 1. 
To each site x G A^,, we will associate a finite measure Hx on C and an element Vx £ B{l-L\j^) 
with the form 

(2.85) Vx= I W{z6x)^x{dz). 

Jc 

Here, for each z £ C, W{zdx) is a Weyl operator as discussed in Section 12.2.21 We require 
that each ^x is even, i.e. invariant under z i— )■ —z, to ensure self-adjointness; namely 
V* = Vx- The anharmonic Hamiltonians we consider are given by 

— ^ 9 9 9 — ^ / \9 — ^ 

= Px+^ Qx + ^j{<lx - Qx+e,) + Z^^x- 

x&Kl i=l a:GAi, 

We denote the dynamics generated by on B{'Hai^) ^ thai is 

(2.87) Tt{A) = e'^^^Ae-''^^ for A G BCHa^) . 
Before we present our Lieb-Robinson bounds, we discuss two examples. 

Example 2.9. Let ^x be supported on [—it, vr) and absolutely continuous with respect to 
Lebesgue measure, i.e. fix{dz) = Vx{z)dz. If Vx is in L^([— vr, vr)), then Vx is proportional to 
an operator of multiplication by the inverse Fourier transform of Vx. Moreover, since the 
support of jjLx is real, Vx corresponds to multiplication by a function depending only on Qx- 

Example 2.10. Let fix have finite support, e.g., take supp{fix) = {z, —z} for some number 
z = a + ip £ C. Then 

(2.88) Vx = W{z6x) + W{-z6x) = 2 cos(Qg^ + Ppx) . 
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We now state our first result. 
Theorem 2.11. Fix L > 1 and take V, Hi, and r/" as defined above. Suppose that 

(2.89) K= sup / \z\'^\iix\{dz) < oo . 

xgAl Jc 

For every fi > and e > 0, there exist positive numbers c and v for which the estimate 

(2.90) \\[rHWif)),W{g)]\\<ce^^+'^-^'^^\'\ \fix)\\giy)\F^ {dix,y)) 

holds for all functions f,g£ £^(Ai) and any t G R. Here 

(2.91) c = C{e,fi) and v = {fi + €)vh{fJ' + e) 

as in Corollary \2.8\ while Cd is the convolution constant for F^ as in \2. 80\) . 

Proof. With L > 1 fixed, we will denote by = t^'^ and = r/' for notational convenience. 
Fix t > and define the function ■ [0, — ?• B{T-L\j^) by setting 

(2.92) ^t{s)=[Ts{rU{W{f))),W{g)\ . 

It is clear that interpolates between the commutator associated with the original har- 
monic dynamics, at s = 0, and that of the perturbed dynamics, rj at s = t. A calculation 
shows that 

(2.93) ^^*(^)=^ E [rs{[V.,W{Tt^sf)]).W{g)] . 
The inner commutator can be expressed as 

[Vx,W{Tt^sf)] = f [W{z6^),W{Tt^sf)]f^Adz) 
Jc 

(2.94) = WiTt-sf)Ct-sAf)- 
where 

(2.95) = ^t-sAf) = ^ Wiz6x) |e^i-[m-./.-^.>l _ i} ^,^(dz) G ^(^aJ . 
Thus satisfies 

^vl/t(.) = i j;^t(.)r, (£,_,;,.(/)) 

(2.96) ^'^"^^ 

+i rs{W{Tt-sf)) [Ts{Ct-sAf))^W{g)] . 

The first term above is norm preserving. In fact, define a unitary evolution Ut{-) by setting 

(2.97) ^Ut{s) = -iY^s {Ct-s-Af)) Ut{s) with Um = t . 

It is easy to see that 

(2.98) ^ {Ms)Ut{s)) =iY^^ {W{Tt-sf)) [r. {CtsM)) , W{g)] Ut{s) , 
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and therefore, 



(2.99) ^titptit) = ^t(O) + i V / Ts {W{Tt^sf)) [ts (Ct-s-Af)) , W{g)] Ut{s) ds . 

J77 -'0 



Estimating in norm, we find that 

'rt {W{f)) , W{g)] II <|| [rO {W{f)) , W{g) 



(2.100) 



+ f [rs{Ct-sAf)),W{g)] 



ds . 



Using Coroharv 12.81 we know that for any /x > and e > 0, 

(2.101) ||[tO(H^(/)),H^(5)]|| < C(e,^)e('^+^)''^('^+^)l*l ^ \fix)\\g{y)\ FAdix,y)) 

Similarly, one can estimate 

\hn[{Tt-sf.z5M < \A\Tt-sf{x)\ 

(2.102) < |z|C(e,M)e(^+^)^'^('^+^)(*-^) ^ \f{x')\FAd{x',x)). 



In this case, the bound 



(2.103) 



[Ts{Ct.sAf))^W{g)] <C(e,^)e('^+^>'^('^+^)(*~^) ^ \f{x')\FAd{x',x)y, 

[Ts{W{z6^)),W{g)] \fi^\{dz) 



X / \z\ 
Ic 



follows from (|2.95|) . Setting c = C{e, fi) and v = (/U + e)t>/j(/i + e), the combination of (|2.100p . 
(12.1011) , and (l2.1(K-;p demonstrate that 



rtiW{f)),W{g)\ <ce^' ^ \f{x)\\g{y)\ {dix,y)) 



x'gAl x&Al 

X / \ A\\[rs{W{z5^)),W{g)] \^l^\{dz)ds. 



(2.104) 



Following an iteration scheme similar to the one in the proof of Theorem 12. 3| one arrives 
at (fZUOD as claimed. □ 

The statement of the Lieb-Robinson bound proven in Theorem 12 . 1 1 1 can be strengthened 
to include a larger class of perturbations. In fact, perturbations involving short range 
interactions can be handled quite similarly. We introduce these perturbations as follows. 

For each subset X C Ki, we associate a finite measure on C-^ and an element 
Vx e B{niy^) of the form 



(2.105) 



Vx= I W{z■5x)^lx{dz), 
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where, for each z £ C''*', the function z ■ bx : — t- C is given by 
(2.106) 



x'ex 



Zx if X G X, 
otherwise. 



We wiU again require that fix is invariant with respect to z i— >• —z, and hence, Vx is 
self-adjoint. In analogy to (I2.86p . we will write 



(2.107) 



V= Vx, 

XcAl 



where the sum is over all subsets of A^. Here, as before, we will let r/" denote the dynamics 
corresponding to + V. 

The main assumption on these multi-site perturbations is as follows. We assume there 
exists a number fii > such that for all < < //i, there is a number /t^ > for which 
given any pair x,y € Al, 



(2.108) 



E 

x,yeX 



kxIkylUxkd^) < K^Ff,{d{x,y)) . 



In this case, the following Lieb-Robinson bound holds. 

Theorem 2.12. Fix L > 1, V , and as above. Assume that \2. 1 08]) holds. Then, for 
any < /i < /ii and e > 0, there exist positive numbers c and v for which the estimate 



(2.109) 



[rtHW{f)),W{g)]\\<ce^^+'--^'^^\'\ J2 \f{x)\\g{y)\F, {d{x,y)) 

x,yGA£, 



holds for all functions f,gG £'^{Al) and any t £ W. Here c, v, and Cd are as in Theo- 
rem \2.n\ 

The proof of this result closely follows that of Theorem 12.111 and so we only comment 
on the differences. 

Proof. For f,ge ^^(Al) and t > 0, define : [0,t] B{nA^) as in I^M^. The derivative 
calculation beginning with (I2.93P proceeds as before. Here 



(2.110) 



Ct-s;Xif) 



W{z-6x)[ 



,ilm[(Tt-s.f,z-Sx)] 



1 1 fixidz 



is also self-adjoint. The norm estimate 

'n {W{f)) , W{g)\ II <|| [rO (Wif)) , Wig) 



(2.111) 



XcA 







[Ts{Ct-sMf)),Wig)] 



ds . 



holds similarly. With (I2.110p . it is easy to see that the integrand in (I2.11ip is bounded by 



(2.112) ce^^'-'^ 1/(^)1 E ^Adi^^^')) [ M 



[Ts{W{z-5x)),W{g)] Midz), 
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the analogue of (|2.103p . for any > and e > 0. Proceeding as before, 
(2.113) 

'rt{W{f)),Wig)\\\<ce''' J2 \f{x)\\g{y)\F^{d{x,y)) 



ft 



xcAl x'ex 



X e^^'~'U \z,,\\\[Ts{W{z-6x)),W{g)] \fix\{dz)ds. 



The estimate claimed in (|2.109p follows by iteration. In fact, the first term in the iteration 
is bounded by 

(2.114) 



o{t-s) 



:€Al XcAl xiex •'^ 

X / kxil(ce'"* ^ ^ \z^^\\g{y)\F^{d{x2,y))^\iix\{dz)ds 

'^^ X2&Xy&AL 

<ct-ce^' \f{x)My)\ E F^{d{x,xi))F^{d{x2.y)) 



x,y£Ai 



x-i,X2£AL 



E / \zxA\zx2\\l^x\{dz) 



x\,X2€X 



<K^ct-ce^' J2 \f{^)My)\ E F^{d{x,xi))F^{d{xuX2))F^{d{x2,y)) 

x,y£AL xi,X2€Al 

< K^Clct ■ ce^* ^|/(x)||5(y)|F^ {d{x, y)) , 



x,y 



where we used that < ^ < /^i in the second inequality above. The higher order iterates 
are treated similarly. □ 



3. Existence of the Dynamics 

The goal of this section is to demonstrate that, in a suitable sense, Lieb-Robinson bounds 
imply the existence of the dynamics in the thermodynamic limit. We prove a general 
statement to this effect in the Section 13.11 for the case of bounded interactions. When 
considering anharmonic systems, more care must be taken in analyzing the thermodynamic 
limit. We discuss recent results in this direction in Section [3. 2[ The analogous problem for 
the classical anharmonic lattice was analyzed in [27j . 

3.1. Bounded Interactions. It is well-known that Lieb-Robinson bounds for quantum 
systems imply the existence of the dynamics in the thermodynamic limit, see e.g. [9l I35j . 
Here we demonstrate that the same basic argument also applies in the general case of 
bounded interactions. Our set-up for this section is similar to that of Section 12. H except 
that now we regard F, which is still equipped with a metric d, as a countable set with 
infinite cardinality. 
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In many examples, our models are defined over F = Z*^ for some d > 1. For locality 
estimates and ultimately a proof of the existence of the dynamics, the underlying lattice 
structure of Z'^ is not a necessary assumption. We express the required regularity of F in 
terms of a non-increasing function F : [0, oo) — t- (0,oo) as mentioned in Section [2.11 

We will say that the set F is regular if there exists a non-increasing function F : [0, oo) — 
(0, oo) which satisfies: 

i) F is uniformly integrable over F, i.e., 

(3.1) ||F|| := supVF(d(x,y)) < oo, 
and 

ii) F satisfies 

(3.2) C := snp ^ ^ W^^^))^ W^. i/)) , 

For finite sets X C F, the Hilbert space Tlx and the local algebra of observables Ax 
are defined as in (|2.1|) and (|2.2p respectively. Recall also that for finite sets X G Y G T, 
Ax C Ay, and we may therefore define the algebra of local observables by the inductive 
limit 

(3.3) = U 

xcr 

where the union is over all finite subsets X gT; see [51 12] for a more detailed discussion. 

Our first result on the existence of the dynamics corresponds to Hamiltonians defined as 
bounded perturbations of local self-adjoint operators. More specifically, we fix a collection 
of on-site local operators H^°'^ = {Hx}xer where each Hx is assumed to be a self-adjoint 
operator over Tlx- In addition, we will consider a general class of bounded perturbations. 
These perturbations are defined in terms of an interaction $, which is a map from the set 
of subsets of F to with the property that for each finite set X C T, ^{X) S Ax and 
^{X)* = ^{X). To prove the existence of the dynamics in the thermodynamics limit, we 
require a growth condition on the set of interactions $ being considered. This condition 
is expressed in terms of a norm analogous to the one introduced in our proof of the Lieb- 
Robinson bounds in Section [2.11 Denote by B(T, F) the set of interactions with 

(3.4) ll^ll := sup } Yl 11^(^)11 < ^■ 

x,yerF{d{x,y)) ^^^^ 

Now, for a fixed sequence of local Hamiltonians H^°'^ = {Hx}xer, as described above, an 
interaction <I> G B{T,F), and a finite subset A C F, we will consider self-adjoint Hamiltoni- 
ans of the form 

(3.5) Ha = H'r + Ht = J2H.+ Yl 

xeA XcA 

acting on T-La (with domain given by ^x€A D{Hx) where D{Hx) C Hx denotes the domain 
of Hx)- As these operators are self-adjoint, they generate a dynamics, or time evolution, 
{r/^}, which is the one parameter group of automorphisms defined by 

T^{A) = e^^^^Ae-'^^^ for any A e Aa. 
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Theorem 3.1 ([37]). Under the conditions stated above, for all t G M and A G Ar, the 
norm limit 

(3.6) lim r,\A) = niA) 

exists. Here the limit is taken along any sequence of finite, non- decreasing, exhaustive sets 
A which tend to T. The limiting dynamics Tt defines a group of *- automorphisms on the 
completion of Ar- In addition, the convergence is uniform for t in a compact set. 

Since the convergence proven in Theorem 13. II above is in norm and the estimates provided 
in Theorem 12.31 are independent of the volume, it is clear that the dynamics defined above 
also satisfies the Lieb-Robinson bound ()2.1ip . 

Proof. Let A C T be a finite set. Consider the unitary propagator 

(3.7) UA{t, s) = e**^^ e-*(*-^)^^ e— 
and its associated interaction-picture evolution defined by 

(3.8) T-f^inti^) = ^a(0. t) ^ ^a(*i 0) for all ^ G • 
Clearly, U\{t, t) = 1 for all t G M, and it is also easy to check that 

t^UA{t,s) = H'f{t)l(A{t,s) and - i-^UA{t, s) =UA{t, s) WJ^\s) 
at as 

with the time-dependent generator 

(3.9) i/f (t) = e^^A^*i/f e-^^^r* = ^ e^^^A^* ^.(Z) e'^^A^* . 

ZcA 

Fix T > and X CT finite. For any A G Ax, we will show that for any non-decreasing, 
exhausting sequence {A„} of F, the sequence {t^^^{A)} is Cauchy in norm, uniformly for 
t G [-T,T]. Since 

an analogous statement then immediately follows for {t^"{A)}. 
Take n < m with X C A„ C A^ and calculate 

(3.10) T,%{A) - T,%,{A) = ")^A„(s, 0)} ds . 

A short calculation shows that 

^Ua^ (0, s) Ua„ {s, t) AUa^ (t, s) Ua,^ {s, 0) 
(3-11) = iUA^{0,s)[{H^^{s) - H^^is)) ,UAAs,t)AUAAt,s)]UAjs,0) 

= ii^A„(0,.)e-^t^ [B{s),Tt, [Ait))] e-<:Z^A„(.,0), 

where 

(3.12) A{t) = e-^*^t^^e**^A„^ = e'^'^'^^e^*^^"^ 
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^A°Vuint^ \ [Tint/ \\ isWP" 
Z<ZKm ZcA„ 



(3.13) 



E 



loc 

Am\An 



ZCAm: 

ZnA„\A„^0 



Combining the results of (j3.10p - p.l3p . we find that 



(3.14) 



< 



ds 



and by the Lieb- Robinson bound Theorem \2.'6\ it is clear that 



(3.15) \[Ttt{A{t)),B{s) 



< 



E 



isH\° 



rf_!'4 (A{t)) ,e A™\A„$(2')e '^'"AmVA" 



ZCAm: 

ZnA™\A„^0 



< 



< 



< 



2PII 

2p|| 

2p||||$| 
C 

< 2||yl||||$|| 



,2||'I>|lC|t-s| 



,2||<I>||C|t-s| 



i) E E 11^(^)11 EE ^(^(^'^)) 



yeAm\A„ ZCAm.: 

yez 



i) E E E 11^(^)11 E^(^(^'^)) 



J/GAm\A„ ZGAm ZCAm: 



,2|l<I>|lC|t-s| 



,2||<I>||C|t-sl 



E E E mx,zmd{z,y)) 

J/eAm\A„ xeX ZGAm 

i) E E^(^(^'y))- 



With the estimate above and the properties of the function F, it is clear that 



(3.16) 



sup 

t€l~T,T] 



as n, m — )• oo. 



This proves the claim. 



□ 



If all local Hamiltonians are bounded, as is the case for quantum spin systems, the 
infinite volume dynamics {r^}, whose existence we proved above, is strongly continuous. If 
the Hx are allowed to be densely defined unbounded self-adjoint operators, we only have 
weak continuity and the dynamics is more naturally defined on a von Neumann algebra. 
This can be done when we have a sufficiently nice invariant state for the model with only 
the on-site Hamiltonians. Suppose, for example, that for each x S F, we have a normalized 
eigenvector (px of Hx- Then, for all A G Aa, for any finite A C T, define 

(3.17) P(^) = ((g)</'x,A(g)0x). 

xeA xeA 
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p can be regarded as a state of the infinite system defined on the norm completion of ^r- 
The GNS Hilbert space Hp of p can be constructed as the closure of Av^^^Y't'x- 
gP 4>x- Then 

\\{rM)-rtM))n < U{A) - tI^-\A)) 
(3.18) ^ ^ 

+ 



(rf ")(^) - r(f")(A)) V|| + II (r,^)(^) - n,{A)) i; 



For sufficiently large A„, the limj_j.jp of middle term vanishes by Stone's theorem. The two 
other terms are handled by (|3.6p . It is clear how to extend the continuity to ip £ Tip. 

3.2. Unbounded interactions. In this section, we will prove the existence of the dynamics 
in the thermodynamic limit for the bounded perturbations of the harmonic Hamiltonian 
we considered in Section 12.2.41 The existence of this limit was considered in a recent work 
[6] where, by modifying the topology, a rigorous analysis of the dynamics corresponding to 
the anharmonic system in the finite volume could be performed in the limit of the volume 
tending to Z'^. Here, as in |37j, we take a different approach. With our method, we regard the 
finite volume anharmonicities as a perturbation of the infinite volume harmonic dynamics. 
We prove that the limiting anharmonic dynamics retains the same weak continuity as the 
infinite volume harmonic dynamics. 

3.2.1. The infinite-volume harmonic dynamics. It is well-known that the harmonic Hamil- 
tonian defines a quasi-free dynamics on the Weyl algebra. We briefly review these notions 
here and refer the interested reader to [9] (see also [37]) for more details. 

In general, the Weyl algebra, or OCR algebra, can be defined over any linear space T> that 
is equipped with a non-degenerate, symplectic bilinear form. For the current presentation, 
it suffices to think of P as a subspace offiZ"^), e.g. £^(Z'^), fiZ"^), or ^^(A) for some finite 
A C Z'^. In this case, the symplectic form is taken to be lm[{f,g)]. 

The Weyl operators over D are defined by associating non-zero elements W{f) to each 
f GV which satisfy 

(3.19) Wif)* = W{-f) for each / G P , 
and 

(3.20) W{f)W{g) = e-^i'^[<-^'f>l/2^(/ + fo^ ah /, 5 € P . 

It is well-known that there is a unique, up to *-isomorphism, C*-algebra generated by 
these Weyl operators with the property that 1^(0) = 1, W{f) is unitary for all f G D, 
and \\W{f) - 1\\ = 2 for all / G P \ {0}, see e.g. Theorem 5.2.8 [9]. We will denote by 
W = W{T>) this algebra, commonly known as the OCR algebra, or Weyl algebra, over T>. 

A quasi- free dynamics on W(P) is a one-parameter group of *-automorphisms Tt of the 
form 

(3.21) rt{W{f)) = W{Ttf), feV 

where Tj : P — )• P is a group of real-linear, symplectic transformations, i.e., 

(3.22) To = 1, Ts+t = Ts o Tt, and, Im [{Ttf, Ttg)] = Im [(/, g)] . 

Since ||VF(/) — VF(i?)|| = 2 whenever f ^ g £ V, such a quasi-free dynamics will not be 
strongly continuous; even in the finite volume. 



22 



B. NACHTERGAELE AND R. SIMS 



To define the infinite volume harmonic dynamics, we must recah the finite volume calcu- 
lations from Section [2.2.2I Let 7 : [— vr, vr)'^ — t- M be defined as in (j2.42p . Take U and V as 
in (j2.6ip with the unitary Fourier transform from to L^([— vr, vr)'^). Setting 



(3.23) 



Tt = {U + V)F-^MtT{U* - V*) 



one can easily verify (|3.22p using the properties of U and V; namely (j2.63p and (j2.64p . 
If, in addition, T> is Tt invariant, then Theorem 5.2.8 of [9j guarantees the existence of a 
unique one parameter group of *-automorphisms on W(I'), which we will denote by r^, that 
satisfies ()3.2ip . This defines the harmonic dynamics on such a yV{T>). 

With calculations similar to those found in p6], one finds that the mapping Tt defined 
above can be expressed as a convolution, analogously to the finite volume calculations. In 
fact. 



(3.24) 
where 

(3.25) 



Ttf = f*[ m 



r(0) 



(-I)- 



H^'\x) 



1 



(27r)'^ 
1 

(2^ 
1 

(2^ 



Im 

Re 
Im 



f J_g^(fc-x-27(fc)i) 
^i{k-x~2i{k)t) 



and we have replaced the Riemann sums from the finite volume with integrals. The following 
result holds. 



Lemma 3.2. Consider the functions defined in 113. 25\) . For w > 0, Ai 
that Ci^^x = (w^ + 4^ . ^j)^^"^ > 0; '^'^^ '^'^y > 0> tfiG bounds 



Aw > 0, but such 



(3.26) 



hI'\x) 



< e 



-A.(|x|-c^,;,max(|,e("/2)+i)|t|) 



,e(''/2)+i)|t|) 
< ;,e^/2g-/<l-h^^,A max(| , e(M/2)+i) 



hold for allteR and x G Z"'. Here \x\ = YlLi 



Given the estimates in Lemma 13. 2| equation (j3.24p , and Young's inequality, Tt can be 
defined as a transformation of P'{'L'^), for p>l. However, the symplectic form limits us to 
consider V = £P{Z'^) with 1 < p < 2. 

Mimicking the arguments from the proof of Theorem 12.61 the above estimates yield the 
following Lieb-Robinson bound for the infinite volume harmonic dynamics rj. 

Theorem 3.3. For any fi > and e > 0, there exist positive numbers c and v for which 
the estimate 



(3.27) 



\\[Tt{W{f)),W{g)]\\<ce^\'\ \f{x)\\g{y)\F,{\ 



X -y\ 
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holds for all functions f,g€z ^ (Z ) and any t G M. Here one may take 
(3.28) c=(l + c^,Ae(^+^)/2 + c-\) sup e-"(l + sy+' 



and 



s>0 



(3.29) v = {fi + e)c^,A max (^-— , e^^+^^/^+i 
mi/ic,,A = (a;2 + 4Ej=iA,)i/2. 

3.2.2. Weak continuity and the anharmonic dynamics. As we indicated in the previous 
section, the harmonic dynamics is not strongly continuous, not even when restricted to a 
finite volume. It is possible, however, to show weak continuity of the harmonic dynamics 
in the GNS-representation of certain states. If /> is a regular, rt-invariant state on W(P), 
then weak continuity follows from proving continuity of the functions 

(3.30) p{W{gi)W{Ttf)W{g2)) for aU 51, 52, / G 2? • 

In [37], we verified these properties for the infinite volume ground state of the harmonic 
Hamiltonian, i.e. the vacuum state for the 6-operators, defined on W(P) by setting 

(3.31) /5(H^(/)) = e-3ll(^*-^*)^ll' , for all /e P. 

Alternatively, one could also prove weak continuity of the harmonic dynamics in a repre- 
sentation corresponding to equilibrium states at positive temperature. In either case, it is 
precisely such a weakly continuous dynamics to which we add our anharmonic perturba- 
tions. 

Using Proposition 5.4.1 from [9], which applies to a general weakly continuous dynamics, 
in fact to a VF*-dynamical system, we define a perturbed dynamics as follows. Fix a finite 
subset A C Z"'. Consider a perturbation of the form Va = ^^eA where, for each x ^ A, 
Vx is as defined in (I2.85P of Section 12.2.41 The arguments below equally well apply to the 
multi-site perturbations, see (|2.105p . considered at the end of Section 12.2.41 however, for 
simplicity, we only state results in the case of on-site perturbations, see |37j for more details. 
Proposition 5.4.1 demonstrates that the Dyson series 

00 „ 

(3.32) rf )(I^(/)) = n{W{f)) + [^*^(^a)' [• • • Km,rt{Wm]] dti ■ ■ ■ dtn 

jQ<ti<t2-<t 

is well-defined. Furthermore, t^^^ is weakly continuous, and there is a consistency in the 
iteratively defined dynamics; j-^^^i^^^^) ^j^^ constructed by perturbing t[^^^ on A2 
given that Ai n A2 = 0. 

As a consequence of ()3.32p . we prove the following Lieb- Robinson bound in [37J. 

Theorem 3.4. Fix a finite set A C Z'^ and let rj:^^ be as defined above. Suppose that 

(3.33) K= sup / \z\'^\fix\{dz) < 00 . 

For every /i > and e > 0, there exist positive numbers c and v for which the estimate 



(3.34) II [rf ) {Wif)),Wig)] || < ce(''+-^'*)l*l ^ |/(x)| \giy)\F, (|x - y\) 

holds for all functions f,g^ £^{Ij'^) and any t G 
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To prove this result, one argues as in the proof of Theorem 12.111 except that the estimates 
from Theorem 13.31 replace those of Corollary 12. 8i The numbers c, v, and Cd, as well as the 
function F^, are exactly as in Theorem 12.111 

We can now state our result on the existence of the anharmonic dynamics. 

Theorem 3.5. Let Tt be the harmonic dynamics defined on W(^^(Z'^)). Take {A„} to be 
any non- decreasing, exhaustive sequence of finite subsets ojU^ . For each x G , let 



(3.35) 

set Va.„ 
(3.36) 



Sx'eA„ '^'^'^ assume that 



sup 

x&d- JC 



\z\\fix\{dz) < oo and 



sup 



\z\'^\fix\{dz) < oo . 



Then, for each f G £^{Z'^) and t G M, the limit 

(3.37) hm "^(W/(/)) 

exists in norm. Moreover, the limiting dynamics is weakly continuous. 

Proof. To show convergence, we estimate \\T^"{W{f)) — t^"^ {W{f))\\, for C A^, by 

This gives 



considering t^" as a perturbation of 



(3.38) r,^"(iy(/)) 
and therefore 
(3.39) 



'iW{f)) + i 



^A„\A,„,rh(W(/)) 



ds , 



T»-(I1'(/))-T^»(»'(/)) 



XGAn\Am 



\t\ 



iW{f)) 



ds . 



Using Theorem 13.41 we find that 



< 



W{z5,),T^^r{W{f)) |/i,|(dz) 



(3.40) 



< c^{v+C.CMt\-s) ^ |/(y)|F^(|y_^|) 



\z\\^ix\{dz) 



Since / G £^(Z'*) and F^ is uniformly integrable, this estimate suffices to prove that the 
sequence is Cauchy. By observation, the proven convergence is uniform on compact t- 
intervals. 

The claimed weak continuity of the limiting dynamics follows by an e/3 argument similar 
to the one provided at the end of Section 13. 1[ □ 

4. The structure of gapped ground states 

4.1. The Exponential Clustering Theorem. The local structure of a relativistic quan- 
tum field theory [16], is provided by the finite speed of light which implies an automatic 
bound for the Lieb-Robinson velocity. This implies decay of correlations in QFT with a gap 
and a unique vacuum [71 SH [T7]. Fredenhagen [15] proved an exponential bound for this 
decay of the form ~ e"'*''^ '^'l , which corresponds to a correlation length of the form ^ < c/7. 
The gap 7 is interpreted as the mass of the lightest particle. In condensed matter physics, 
the same relation between the spectral gap and the correlation length is widely assumed. 
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The role of the speed of hght is played by a propagation speed relevant for the system at 
hand, such as a speed of sound. A strict mathematical relationship, however, only holds in 
one direction: a unique ground state with a spectral gap implies exponential decay of spatial 
correlations under quite general conditions, which in particular imply a finite bound on the 
speed of propagation known as a Lieb-Robinson bound. This was proved only relatively 
recently [MlESj, using an idea of Hastings [21J. 

As a consequence of subsequent improvements of the prefactor of the Lieb-Robinson (see 
|40j ) ■ we now also have better constants in the Exponential Clustering Theorem than in 
the first results. In particular, for observables with large support it is significant that the 
prefactor is only proportional to the smallest of the surface areas of the supports of the two 
observables. E.g., this is important in certain applications (see, e.g., [211 [33]). 

Theorem 4.1 ([iQj). Let ^ be an interaction with ||$||a < c« for some a > 0. Suppose H 
has a spectral gap 7 > above a unique ground state {•). 

Then, there exists /x > and a constant c = c(-F, 7) such that for all A G Ax, B G Ay, 

\{AB) - {A) {B)\ < c\\A\\ \\B\\mm{d^X,d^Y)e-^"^^^'^\ 

One can take 

a7 



Using the Lieb-Robinson bounds for oscillator lattices one can also prove an exponential 
clustering theorem for these systems. 

Theorem 4.2 ([36]). Let H be the anharmonic lattice Hamiltonian with A > satisfying 
the conditions of Case (ii), and suppose H has a unique ground state and a spectral gap 
7 > above it. Denote by {■) the expectation in the ground state. Then, for all functions 
f and g with finite supports X and Y in the lattice, we have the following estimate: 

\{WU)W{g))-{W{f)){W{g))\ 

<C||/|Ulblloomin(|X|,|y|)e~'^(^'^)/« 

where E, = (4af -|- 7)7(07) and, if we assume d{X,Y) > ^, C is a constant depending only 
on the dimension v. 

The central argument in the proof of these theorems is the same. Here, we only provide 
a sketch and refer to [23l [38l HO] and [36] for the details. 

Suppose H > with unique ground state ii, HVt = 0, with a gap 7 > above 0. Let 
A^Ax and B ^ Ay, d{X, Y) > 0, a,C,v > 0, such that 

\\[Tt{A),B]\\ < C\\A\\ ||5||e-<^Wx.n-^|i|) . 

We can assume {^},AQ) = {i},Bi}) = 0. We want to show that there is a ^ < 00, indepen- 
dent of X,Y,A,B, s.t. 

\{n,ABn) \ < Ce-'^(^'^)/^ 

For 2; € C, Imz > 0, define 

POO 

f{z) = {n,AT,{B)n) = / e''^d{A*n,PEBn) . 

J-y 

For T > 6 > 0, and Tt the upper semicircle of radius T centered at 0: 

2in Jr-r z - lb 
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Then 



\{n,ABn)\ < limsup — 

f)4.o,Ttoo 27r 



T 



fit) 
t - ib' 



dt 



Next, introduce a Gaussian cut-off and remember f{t): 



\{n,ABi})\ < limsup 



1 

2^ 



-at- 



{n,ATt{B)n) 

t - ib 



dt 



assuming 7 > 2ab. For a{d{X,Y)/v)'^ ^ 1, the Lieb-Robinson bounds lets us commute 
Tt{B) with A in this estimate. Using the spectral representation of r^, we get 



\{n,ABn)\ < limsup — 
HO,Ttoo 27r 



T -iEt -ai^ 

dt — d{B*n, PeAVL) 

t — ib 



+ err. 



-7V{4a) 



Optimizing a gives the bounds 



The t— integral can be uniformly bounded by e" 
stated in Theorems 14.11 and 14.21 

The condition that the ground state be unique can be relaxed. E.g., one gets the same 
result for each gapped ground states of infinite systems with several disjoint ground states. 
One can also derive exponential decay in the average of a set of low-energy states separate 
by a gap from the rest of the spectrum, if the number of states in the set does not grow too 
fast with increasing system size. Another straightforward extension covers models of lattice 
fermions [23]. 

The Exponential Clustering Theorem says that a non-vanishing gap 7 implies a finite cor- 
relation length ^. But one can say more about the structure of the ground state. Motivated 
by the goal of devising better algorithms to compute ground states and questions related to 
quantum computation, a number of further results have been derived. The best known is 
Hastings' proof of the Area Law for the entanglement entropy in one dimension [22j, which 
used an approximate factorization lemma of the ground state density matrices. Before we 
discuss this result and a generalization of it, we make a small detour to Valence Bond Solid 
(VBS) models and Matrix Product States (MPS). VBS models were first introduced by 
Affleck, Kennedy, Lieb, and Tasaki [H [3] . MPS are a special case of Finitely Correlated 
States [Ti] . 

The first and best known VBS model is the AKLT model named with the initials of its 
inventors. This model, itself motivated by Haldane's work |18tll9j. led to a dramatic change 
in our outlook on quantum spin chains and the ground states of quantum spin Hamiltonians 
in general. Before the AKLT model, practically all our understanding of the ground states 
of quantum spin systems stemmed directly from exact solutions of special models, primarily 
Bethe-Ansatz solvable models. The Bethe-Ansatz solutions are tremendously important in 
their own right but they had seriously biased our thinking about more general models. 
The AKLT model and subsequent generalizations changed that and led to the much better 
understanding of generic behaviors of quantum spin systems that we now have. So, a small 
excursion to the AKLT model is certainly justified. 

4.2. The AKLT modeL The AKLT model is a spin-1 chain with the following isotropic 
nearest-neighbor Hamiltonian: 



(4.1) 
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acting on y-ia^b] = (^c3-jiX)(6-a+i)^ where is the vector of spin-1 matrices acting on the rrth 
factor. A straightforward computation based on the representation theory of SU(2) shows 
that 

- + -Sx ■ Sx+l + -{Sx • Sx+l)"^ = Px^l+i 

(2) 

where P^. is the orthogonal projection onto the spin-2 subspace of two spin I's at x and 

X + 1. Therefore, -H'i^^ > 0. As we will see in a moment, dim ker f,] = 4, for all 
a < b. Hence, the ground state energy of the model vanishes for all finite chains. 

The AKLT chain has the three properties that characterize the so-called Haldane phase: 

• It has a unique ground state for the infinite chain. In particular, for L > 1, pick 
ipL ^ ker , with HV'lII = 1- Then, for all finite X and A £ Ax, one has a 
limiting expectation value 

u{A) = lim {ijL,Ai>L) 

L—^oo 

which is independent of the chosen sequence. It follows that u is a translation and 
SU(2) invariant state of the quasi-local algebra of observables of the infinite chain. 

• The unique ground state w has a finite correlation length: there exists ^ > 0, C > 0, 
s.t., for ah A £ Ax,B £ Ay 

\u{AB) - io{A)uj{B)\ < C7p||||5||e-'^(^'^)/«. 

In fact, the bound holds with e~^/^ = 1/3 and is optimal. 

• The AKLT chain has a spectral gap above the ground state: there exists 7 > 0, 
such that for all b > a, the gap of H^^^b], which equals the smallest strictly positive 
eigenvalue Ei, satisfies -Ei > 7. For the infinite chain this is expressed by 

uj{A*HxA) > juj{A*A). 

for all X and all A G Ax, with 

{x,x + l} 

Using the Density Matrix Renormalization Group (DMRG) [51], one can compute 
7 numerically to virtually any desired accuracy. E.g., Huse and White found 7 ~ 
.4097... [50]. It was quickly understood that the DMRG can be understood as 
a variational approximation using Matrix Product States (MPS). Since MPS are 
dense in the set of all states [13], the error of this approximation can, in principle, 
be made arbitrarily small. See [5l| for a detailed discussion of the DMRG and 
|451 ITT] for a recent reviews. 

Haldane predicted these properties for the integer-spin Heisenberg antiferromagnetic 
chains. A proof of the existence of non-vanishing spectral gap, or even of the (slightly) 
weaker property of exponetial decay of correlations in the ground state of the (standard) 
Heisenberg quantum spin chains has so far proved elusive, although some interesting con- 
ditional statements have been obtained [Hll]. 

The AKLT chain was the first proven example of the existence of the Haldane phase. 
This is important, but the impact of the explicit construction of the exact ground state of 
the AKLT Hamiltonian has gone a great distance beyond that example. It led to analytic 
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and numerical techniques to compute and approximate the complex entangled states that 
occur in many condensed matter systems (see, e.g., |48i I49j ). 

4.2.1. The AKLT state and its properties. Recall the Clebsch-Gordan series for the decom- 
position of the tensor product of two irreducible representations of SU(2): 

Let (/) S be the singlet state given by 

0=i=(|ti)-|lt», 

and let TV : — )• (8) be the isometry implementing the embedding corresponding to 

any observable of the spin-1 system at a single site, A G M3, 
WAW* is its embedding in M2 (8) M2. Then, for every n > 1, and any |a) , G C^, define 

(n) 

the vector V'^^/i G '^[i.n] by 

(4.2) = {W* ■ ■ ■ <^W*){\a) ^ (f)<^ ■ ■ ■ (P0\I3)). 

Since W* intertwines SU(2) representations, so does W* ■ ■ ■ W*. In particular, W* <^ 

(n) 

• • • (8) W* leaves the total spin of any vector unchanged. Since ip^^ is the image of a vector 

spin does not exceed 1. It follows immediately that V^lg is a 
ground state of H[i^n]j because -ff[i,n] is a sum of projections on the spin-2 states of a pair 
of neighboring spins: 

Pi,%,(.W*0W*){\a)(S>cP(S>\/3)) = O 

It is not hard to show that the vectors of the form ip'^^ in fact span ker H[i^n]^ i-^-, all ground 
states of -?^[i,n] are of this form. 

To show the uniqueness of the thermodynamic limit and the finiteness of the correlation 
length, we consider the structure of the expectation of an arbitrary observable: 

MAi<^---(S)A^)= ^. 

Careful inspection reveals that we can write this formula in the following form 

(4.3) a;„(>li ® • • • ® ^„) = C7„T:P„Ea, o Ea^ o • • • Ea^P^) 
where, for A e M3 and B e M2, E,a{B) G M2 is defined as 

Ea{B) = V*A(S)BV 

with y : C2 ^ C3 C2 defined by 

V\a) =c(iy*® ll2)(|a) 0). 

and c and Cn are normalization constants. It follows from the properties of the singlet 
vector (p and the intertwining operator W* , that V is also an intertwiner. By choosing the 
constant c we can make V the up to a phase unique isometry corresponding to the inclusion 
15(1/2) ^ ^(1) £)(i/2)^ -^i^j^ ^Yiis choice it is clear that 

Ei3(ll2) = I2. 
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The normalization constant C„ in (|4.3p is then simply equal to 1. One can further check 
by a simple computation that 

which is equivalent to the statement that the linear map Ejg is diagonal in the basis of AI2 
consisting of II2 (with eigenvalue 1) and the three spin- 1/2 matrices (each with eigenvalue 
—1/3. The kth powers of the Ejg are therefore given by 

(4.4) {^i,)HB) = \{TtB)1, + {-l^iB - ^TtB), 

for all B G M2. From this property it follows immediately that the thermodynamic limit of 
the formula (j4.3p exists and is independent of the choice of Pa and P^: 

(1 • • • 1 ^Ai (g) • • • ^„ (g) 1 • • • (g) 1) 

n; nr 

(4.5) =^TvEa,oEa,o...EaM 

It is also clear that the convergence is exponentially fast and, by the same consideration, 
that the correlations in this state decay as (i/3)distance_ 

The third essential property of the AKLT model is the non-vanishing spectral gap. In 
view of our later discussion of the Area Law and approximate factorization property of 
gapped ground states, it useful to present to underlying structure of the AKLT ground 
state in a bit more detail. 

Let u denote the ground state of the infinite AKLT chain defined by (j4.5|) . Let P[a,b] be 
the density matrix describing the restriction of u to ^[a,fe], i-e., for all Aa, . . . ,Ah £ M3, 

ujniAa ^■■■(S)Ab) = Trp[^^^Aa ® • • • ® A. 

Then, the rank of p[a,b] is 4 (equal to the nullity of ffj^^ ). Let G^a,b] be the orthogonal 
projection onto the range of pj^ ^ . Then, again using (j4.4p . one can show that for £ > 0, a > 
i + l 

(4.6) \\G[a-e,a+e+i] [G[i,a] ® G[a+i,L]] - G[i^L]] II < Ce-^/^ 

This property allows one to prove a uniform lower bound for the gap \14:\ [Ml B6] . In brief: 

7 > 1(1 - ce-^/€) X (gap of fAKLT)^ 

4.3. The Area Law for the entanglement entropy. The AKLT state u satisfy an "area 
bound" on the entropy of its local restrictions. In general, this means that for X C A and 
px G •^x is the density matrix describing the restriction of the state to Ax, then 

(4.7) S{px) = -Trp^log px<C\dX\ 

If X is an interval, for the AKLT state we have S{px) = log 4. This is a trivial consequence 
of the fact that 4 is the rank of px- The Area Law Conjecture: (|4.7p holds in general for 
gapped ground states of arbitrary quantum spin systems with bounded spins and bounded 
finite-range interactions. The conjecture has been proved by Hastings in the case of one- 
dimensional systems [22J. The evidence for the higher dimensional case comes primarily 
from the existence of a large class of models (VBS models) with AKLT-like ground states 
in arbitrary dimensions and their ground state are sometimes called PEPS (Products of 
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Entangled Pairs) [39]. In one dimension there is a density result stating that VBS states 
with a similar structure as the AKLT state are weakly dense in the set of all pure translation 
invariant states [13]. It is also know that each such state is the unique ground state of 
a finite-range Hamiltonian with a non- vanishing spectral gap [H]. If we assume that a 
similar genericity holds for the higher dimensional VBS states, the evidence for the Area 
Law Conjecture is quite strong. The rank of the local density matrices in a VBS state is 
bounded by the dimension of the space of boundary vectors. In one dimension these are 
the vectors |a) and that appear in (|4.2p . This dimension is of the form d^^'^^, leading 
immediately to a bound of the form (j4.7p . 

The theory of higher-dimensional VBS models is still in progress. In the next section we 
present a result that is consistent with the assumption that the unique gapped ground states 
of finite-range Hamiltonians in higher dimensions may indeed be well approximated by VBS 
states. Specifically, we will see that gapped ground states in general have an approximate 
product structure similar to (j4.6p . 

4.4. An approximation theorem for gapped ground states. We will consider a sys- 
tem of the following type: Let A be a finite subset of Z"^. At each x G A, we have a 
finite-dimensional Hilbert space of dimension Ux- Let 

Hv= Yl ^ix,y), 

{x,y}cA,\x-y\=l 

with ||<I>(x,y)|| < J. Suppose Hy has a unique ground state and denote by Pq the corre- 
sponding projection, and let 7 > be the gap above the ground state energy. 
For a set vl C A, the boundary of A, denoted by OA, is 

dA = {x €z A\ there exists y £ A\A, with Ix — y| = 1}. 

and for i > 1 define 

B{e) = {x£A \ d{x,dA) < £}. 

The following generalizes to arbitrary dimensions a one-dimensional result by Hastings 
[22]. 

Theorem 4.3 ([20]). There exists ^ > (given explicitly in terms of d, J, and 'j), such 
that for any sufficiently large m > 0, and any A C A, there exist two orthogonal projections 
Pa S Aa, and Pa\a ^ '^A\A; o.nd an operator Pb G ^B{m) ^^^^ II-P-bII ^ 1; such that 

\\Pb{Pa ® Pa\a) -Po\\< CiOm'e-"'/^ 
where C(^) is an explicit polynomial in ^. 

The proof of this theorem uses several ideas of [22]. Below, we only outline three main 
steps and refer to |20] for the details. 

(1) The first step is the bring the Hamiltonian in a form similar to the Hamiltonian of 
the AKLT model in the sense that we split the Hamiltonian into terms (three in this case) 
which are each individually minimized by the ground state, up to some error we can make 
arbitrarily small. The three terms correspond to the set A and its complement A \ ^, and 
a boundary of thickness i to describe the interaction between A and A \ ^4. By taking i 
sufficiently large the error can made arbitrarily small. 

Without loss of generality we can assume that the ground state energy of H\ vanishes: 
H^ipo = 0. We aim at a decomposition of H\, for each sufficiently large £, into three terms: 

Ha = Ka + Kb[i) + Kx\A ! 
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Figure 1. The set A and its fattened boundary B{£) used in the statement 
of Theorem 14.31 

with the following two properties for each Kx, X = A, B(i),A \ A: 

(i) suppKx C X; 

(ii) ll-fTx^oll < e~^^, for each X and for some c > 0. 

Note that we only assumed H^ipQ = 0, and no special properties of the interaction terms 
We start from 

(4.8) Ha = Hj + Hb + He, 

where 

/ = I{e) = {x G A I for all y G dA, d{x, y) > £} 
E = E{£) = {x eV\A\ioT ally e dA, d{x, y) > £} . 
The sets I{£) and E{£) are the interior and exterior of A. B{£) is boundary of thickness 2£: 

B{£) = {x G A I d{x,dA) < £} . 
Note that A is the disjoint union of /, B, and E. Now define 
Hi= ^{X), Hb=Y.HX), He= 

XCA: XCA: XCA: 

For £ > 1, there are no repeated terms and ()4.8p holds. However, there is no guarantee 
that ll-ffjc^oll will be small. In general, this will not be the case but we can arrange it so 
that each term has expectation in ^o- What is needed is a bit of 'smoothing' of the terms 
using the dynamics: for X G {/, B, define 

{Hx)a = \[^ r Tt{Hx)e--''dt, 
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for a > 0. Since the full Hamiltonian is invariant under the dynamics it generates, we still 
have 

Hk = {Hi)a + {HB)a + {HE)a , 

But now we can show that HffxV'oll is small for a small. Unfortunately, the support of 
{Hx)a is no longer X. The easiest way to correct this is by redefining them with a suitably 
restricted dynamics as follows: 

K^f = e^'^'^Hje-^'''^ e""*' dt , 

and similarly define k'^^^j^ using Hx\a, and K^^^ using -ffs(2^)- 

A good choice for a is av^ /{2i), where a and v are the constants appearing in the Lieb- 
Robinson bounds for the model under consideration. With this choice one can show that 
all errors are bounded by 

e(^) = C{d,a,v)j'^\dA\l'^-^l'^e-^'^ 

with 

^ = 2 max(a~^, af ^/7^) 
To summarize, in step (1) we obtained an approximate decomposition 

\\H^-{Kf + KP + K'^^:^^)\\<e{l) 

with the desired property 

||i^?Vo|| <e(^) 

for X = A,B,K\A. 

(2) Next, we define the projections Pa and P\\a as the spectral projections of K^^'^ and 
-^A\A projecting onto their eigenvectors with the eigenvalues < ^/eJJ). This gives 

\\{t-PA)M\ < -^WK^A^M < vW) 

and similary for P\\a- Since the projections commute we have the identity 

2(1 - PaPa\a) = (1 - ^a)(1 + Pa\a) + (1 - ^a\a)(]1 + Pa) , 
from which we obtain 



(4.9) ||Po - PoPaPa\a)\\ = \\Po{T^ - PaPa\a)\\ < ■ 

(3) As the final step, we need to replace the ground state projection Pq which multiplies 
PaPa\a in the LHS of (14. 9p . by an operator supported in the boundary set B{m) for a 
suitable m. 

We start from the observation that for a self-adjoint operator with a gap, such as H\, 
the ground state projection Pq can be approximated by Pa defined by 



Pa = J- I e''"^-e-"' dt. 
vr 



oo 



If the gap is 7, and with our choice of a, we have 

||P^_Pg|| <e-7V(4«) <e-£/? 

We modify this formula for Pa in two ways: 
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(i) we replace e'*^"^ by 



This leads to an operator Pb defined by 



(ii) Then we approximate the result by an operator supported in B(3i) to obtain our 
final results for Pg which appears in the statement of the theorem: 

Pb = TiUf,\B{^,)PB ■ 

With these definitions it is straightforward to show that both ||Po-P4-fA\A ~ PbPaPk\a\\ 
and \\Pb — Pb\\ are small. This concludes the outline of the proof of Theorem 14.31 



5. Conclusions and Further Developments 

In these lecture notes we have reviewed the derivation of Lieb-Robinson bounds for a 
considerable variety of systems, including many of the well-known models frequently used 
in condensed matter physics. It is important to continue to expand the class of systems 
for which Lieb-Robinson bounds can be proved. Their relevance keeps growing as new 
applications continue to be found. 

An application we have not discussed in these notes is the higher-dimensional version 
of the Lieb-Schultz-Mattis Theorem. The classical Lieb-Schultz-Mattis Theorem |28j is for 
spin-1/2 spin chains and states that if the ground state is unique, then the gap above it 
must vanish at least as C/L, where C is a constant and L is the length of the chain. Later, 
Affleck and Lieb generalized the result to other one and quasi-onedimensional models [1]. 
In particular, their result applies to those chains of even length with spins having arbitrary 
half-integer magnitude. But it took more than forty years for someone to make real progress 
on a higher-dimensional analogue of the Lieb-Schultz-Mattis Theorem. In 2004 Hastings 
found a novel approach using Lieb-Robinson bounds directly and indirectly (through the 
Exponential Clustering Theorem), that allowed to obtain a Lieb-Schultz-Mattis theorem in 
arbitrary dimension [21j . The result applies to a wide class of Hamiltonians, which includes 
the half-integer spin antiferromagnetic Heisenberg model on with suitable boundary 
conditions and states that if the ground state is non-degenerate the gap of the system of 
linear size L, 7^, must satisfy: 

(5^1) < C'^^. 

The detailed conditions of the theorem and a rigorous proof are given in [39] . An overview 
can be found in [30] . 

A new application of Lieb-Robinson bounds and their consequences, which recently ap- 
peared on the arXiv, is concerned with the Quantum Hall Effect [24j . This work is concerned 
with system defined on a two-dimensional lattice with torus geometry, with interactions that 
are uniformly bounded and of finite range, and which preserve charge. The authors of [24] 
prove that if such a system has a unique ground state with a non- vanishing spectral gap, 
7, above it, its Hall conductance, (Tx,y, as defined by the Kubo formula, will show sharp 
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quantization. More precisely, for a system of linear size L, an estimate of the following type 
is obtained. There is an integer n, and constants C and c > 0, such that 

(5.2) 



x,y 



n- 



h 



< CL^e 



-C72L2/5/(logL)6 



An even more recent application is the stability of the Toric Code model [26] under small 
perturbations of the interaction [10]. This result significantly enhances the plausibility of 
implementing quantum computation using topologically ordered ground states. Again Lieb- 
Robinson bounds and its corollaries play a crucial role in turning "adiabatic continuation" , 
a tool pioneered by Hastings [2H 125] . into a practical tool for the proof of this result. 
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